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ABSTRACT 

Electromagnetic (EM) counterparts to supermassive black hole binary mergers observed by LISA can be lo- 
calized to within the field of view of astronomical instruments (~ lOdeg^) hours to weeks prior to coalescence. 
The temporal coincidence of any prompt EM counterpart with a gravitationally-timed merger may offer the best 
chance of identifying a unique host galaxy. We discuss the challenges posed by searches for such prompt EM 
counterparts and propose novel observational strategies to address them. In particular, we discuss the size and 
shape evolution of the LISA localization error ellipses on the sky, and quantify the corresponding requirements 
for dedicated EM surveys of the area prior to coalescence. The basic requirements of a wide field of view and 
fast detectors are similar to those for searches being planned for distant cosmological supernovae. A triggered 
EM counterpart search campaign will require monitoring a several-square degree area. It could aim for vari- 
ability at the 24-27 mag level in optical bands, for example, which corresponds to 1-10% of the Eddington 
luminosity of the prime LISA sources of ~ (lO^-lO^)Mo BHs at z = 1 -2, on time-scales of minutes to hours, 
the orbital time-scale of the binary in the last 2-A weeks of coalescence. A cross-correlation of the period of 
any variable EM signal with the quasi-periodic gravitational waveform over 10-1000 cycles may aid the de- 
tection. Alternatively, EM searches can detect a transient signal accompanying the coalescence. The triggered 
searches will be ambitious, but if they successfully identify a unique prompt electromagnetic counterpart, they 
will enable new fundamental tests of gravitational physics. We highlight the measurement of differences in the 
arrival times of photons and gravitons from the same cosmological source as a valuable independent test of the 
massive character of gravity, and of possible violations of Lorentz invariance in the gravity sector 
Subject headings: 



1. INTRODUCTION 

The detection by the future Laser Interfewmetric Space 
Antenna {LISA) of gravitational waves (GW) emitted during 
the coalescence of supermassive black holes (SMBHs) in the 
mass range ^ (1O"^-1O^)M0/(1 +z) will constitute a milestone 
for fundamental physics and astrophysics. High angular res- 
olution observations of the nuclei of several dozen galaxies 
over the past several years have revealed that supermassive 
black holes (SMBH s) are ubiquitous in local gala xies (see, 
e.g., the review by iKormendy & Richstone I Il995h . Merg- 
ers are expected to play an important role in the formation 
and evolution of SMBHs on cosmological time-scales. The 
number of events that will be detectable by LISA is uncer- 
tain by orders of magnitude, depending on model assump- 
tions (Haehnelt 1994; Men ou et al. 1 1200 It IWl onteri et al 



(e.g. IVecchioll2004t iLang & Hughes II2006I) . In addition to 
providing crucial new information for hierarchical structure 
formation scenarios and black hole astrophysics, the obser- 
vation of these GW sources would offer new tests of general 
relativity and cosmology (Hug hes & Menoul2005l:lBerti et al.l 
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2003; Micicetal. 2007; Sesana. Volonteri. & HaardU i2007l 
Islam, Taylor & Silk 2004j). However, for many events in the 
above mass and redshift range, LISA would be able to mea- 
sure the masses and spin vectors of the SMBHs, their or- 
bital parameters, and their luminosity distance, to a precision 
unprecedented in any other type of astronomical observation 
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relativity and cosmo logy (Hugl 
i2005aiialArun et al.ll2006a.b) . 

While the GW signatures themselves are a rich source 
of information, the secure identification of the LISA source 
in electromagnetic (EM) wavebands would enable a host of 
new applications. Traditional astronomical tools could be 
used to study the nature of these sources in detail, clarifying 
how the SMB Hs grow, evolve, an d affect their galaxy envi- 
ronments (e.g. ' Kocsis et al.ll2006h . The simultaneous study 
of photons and gravitons from a single source would also 
likely p robe fundamental aspects of gravitational physics (e.g. 
iDeffayet & Menou 2007) . 

The single most important unresolved issue is whether 
counterparts to LISA sources will, in fact, be found in elec- 
tromagnetic (EM) wavebands. The final angular localization 
by LISA itself is relatively poor (typically a few x 0. 1 square 
degrees at merger), and early studies on counterparts pointed 
out that this large solid angle will contain many thousands 
of galaxies, suggesting that identifying the galactic host of 
the coale scing SMBHs will be difficult, or impossible (e.g. 
IVecchio 1 12004ft . However, recent work has reached much 
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more optimistic conclusions. Using LISA'S localization of the 
source within a 3-D error volume (rather than 2-d sky posi- 
tion) reduces the n umber of ca ndidate galaxi es by ~ 2 - 3 or - 
ders of mag nitude dHolz & Hu ghes 2005; Ko csis et al.ll2006l) . 
Furthermore, if the coalescing SMBHs themselves produce 
bright emission (such as bright quasars), or are associated 
with some other, similarly rare subset of all galaxies (such 
as ultra-luminous infrared galaxies [ULIRGs], whose large 
infrared luminosity may arise from the galaxy merger process 
that also causes the SMBH binary coalescence), then the iden- 
tification of a unique counterp art may become feasible for a 
typical LISA source CiKocsis et al. 2006). 

Another possibility, raised by Kocsis et al.l (l2007l hereafter 
Paper I), is to monitor the sky for EM counterparts in real 
time, even as the SMBH inspiral proceeds. Various strategies 
to carry out this monitoring effort are the subject of the present 
follow-up paper. If the exact nature of counterparts or host 
galaxies were securely known from ab-initio models, then 
one could hope to identify them post-merger, e.g. based on 
the peculiar photometric or spectroscopic signatures. How- 
ever, such ab-initio predictions are difficult. Searching for 
a variable or transient EM signal produced during the timed 
GW-emitting phase could be the most efficient way to iden- 
tify a secure EM counterpart. 

The presence of gas is believed to promote the rapid co- 
alescence of a SMBH binary (see, e.g ., Escalaetal. 2004, 
|2005[ but s ee Milosavlievic & Merritt I i2003 and the review 
by iGualandris & Merritt ,20071 for the possibilitv of coales- 
cence without gas, through stellar dynamical processes). If 
one or both SMBHs continue(s) to accrete gas in the GW- 
emitting stage, for i nstance, the coalescing bin ar y may shine 
as a b right quasar (iKocsis et al.il2006t lArmitage & NataraianI 
l2002h with potential variability on a timescale of ~ hours 
to days, owing to the inspiraling orbital motion of the bi- 
nary. This differs from post-merger counterpart scenarios. 
For example, the merger of two nearly equal-mass SMBHs 
embedded in a thin accretion disk with a central gap may 
be expected to produce an X-ray afterglow, delayed by a 
disk viscous time relative to the merger, i.e. possibly > 
years later (iMilosavljevic & Phinnev i r2005'). Importantly for 
prompt counterpart scenarios, the power emitted in GWs dur- 
ing the final hours of the coalescence is very large, exceeding 
10~^c^/G ~ 10^^ erg s~'. As a result, any coupling of even a 
small fraction of this unprecedentedly large power to the sur- 
rounding gas could produce a luminous counterpart, possibly 
in the form of a variable EM source. 

Motivated by these considerations, in Paper I we devel- 
oped a formalism (which we call the Harmonic Mode Decom- 
position, or HMD method) to efficiently compute the time- 
dependent 3-D localization information contained in the LISA 
data stream. We presented an in-depth study of the potential 
for pre-merger localizations over the large parameter space of 
potential sources. In particular, we showed that using the GW 
inspiral signal accumulated up to some look-back time, ff, 
preceding the final coalescence, one akeady has a good par- 
tial knowledge of the sky position and radial distance to the 
GW source. Indeed, since the sky position is deduced pri- 
marily from the detector's motion around the Sun, we found 
that angular positioning uncertainties are such that often a 
source can be localized to within a few square degrees sev- 
eral weeks prior to merge r. These conclus ions have been con- 
firmed in a recent study ( Lang & Hughes 2007) that included 
the additional effects of spin precession and parameter cross- 



correlations. These effects, which are ignored in our work, 
impact localization errors only several days or less prior to 
merger (resulting in significant improvements at late times). 

The main proposal put forward here is to search for an EM 
counterpart by monitoring LISA source candidates in the ^ 2- 
3 weeks preceding the merger We emphasize that in this con- 
text, the relevant criterion is not the comparison between the 
number density of potential counterparts and the LISA error 
box (^ 0.1 deg^ in a typical case), but the much less stringent 
requirement that the LISA uncertainty be comparable to the 
field of view (FOV) of a wide-field survey instrument (say, 
~ lOdeg^, such as planned for the Large Synoptic Survey 
Telescope, LSST)'. If this requirement is met, it would be 
possible to use such an instrument to monitor LISA coun- 
terpart candidates in real time, during the final stages of the 
BH merger The goal of the present paper is (i) to discuss, 
quantitatively, the feasibility of such a prompt EM counter- 
part search, and (ii) to identify the issues that will likely drive 
the observational strategies. In particular, we derive source 
position uncertainties as a function of time using the results 
of Paper I, and investigate survey strategies to best utilize this 
localization information. 

The remainder of this paper is organized as follows. In 
§ 12] we summarize and expand on the results of Paper I and 
present time-dependent LISA sky position errors. In § |3] we 
discuss the possibility of identifying the EM counterpart of 
the LISA source by monitoring the merger with a wide field 
astronomical instrument. In § H] we discuss the additional 
science that would be allowed by a real-time EM counterpart 
identification (as opposed to identifying a counterpart post- 
merger). In § |5] we summarize the implications of our find- 
ings and conclude. 

2. TIME-DEPENDENT LISA SKY POSITION ERRORS 

The first and most fundamental question in searching for 
prompt EM counterparts is the accuracy to which the LISA 
source can be localized at various look-back times prior to the 
coalescence. In this section, we present these time-dependent 
uncertainties. We use the Harmonic Mode Decomposition 
method described in Paper I. The reader is referred to that pa- 
per for details about the method; here, for completeness, we 
summarize the main assumptions before describing relevant 
results. 

2.1. Assumptions 

1. In general, an SMBH inspiral is described by a total 
of 17 parameters. These include 2 redshifted mass pa- 
rameters,^ (A^2, ?/), 6 parameters related to the BH spin 
vectors, Pspin, the orbital eccentricity, e, the source lu- 
minosity distance, d^, 2 angles locating the source on 
the sky, {6^, (p^), 2 angles that describe the relative ori- 
entation of the binary orbit, {Onl, (^nl), a reference time 
and a reference phase at merger, i.e. at the innermost 
stable circular orbit (ISCO), fisco, ^^'iscOj and, finally, 
the instantaneous orbital phase, (/)oib (which specifies 
the time-evolution of the binary). We denote the look- 
back time before merger by ff. Throughout the paper, 

' Since the redshift of the source can also be inferred from the LISA dis- 
tance measurement, it helps if at least photometr ic reds hifts can be measured 
for all candidates; we will discuss this issue in § l3.2.2l below. 

^ For component masses mi and m2, the total mass is M = m\ +1112, the 
reduced mass is /i = m\m2/M, the symmetric mass ratio is 77 = ^l|M and 
the chirp mass is defined as = Mrf"!^ jMisner. Thome. & WheeleJ1973ll . 
with redshifted masses written as = (1 +z)M- 
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all time parameters will be quoted in redshifted units 
(i.e. as measured by an observer on Earth). 

2. The signal for a GW inspiral is determined by the above 
set of parameters and two additional angles describing 
the orientation of the LISA constellation, (S, $). Note 
that these angles precess with time; indeed, this time- 
dependence carries all of the sky localization informa- 
tion. Since the period is precisely equal to LISA'S or- 
bital period of one year, we need to specify only some 
reference values (e.g. chosen at the time of the merger, 
Sisco and $isco)- 

3. We use the circular, restricted post- Newtonian (PN) ap- 
proximation for the GW waveform (ICutler & Flanagan! 
[1994) keeping only the quadrupole GW harmonic for 
the GW amplitude. Higher order GW harmonics are 
negligible if the quadrupole harmonic is within the 
LISA-frequency band for at least a few months, but 
can improve uncertainties for large mass binaries with 
M;(77/0.25 )^/^ > 2 X IO^Mq an d much shorter observa- 
tion times (HeliTng s & Moorell2003l; I Arun et al.|[2007al: 
Urias & Sintes 20071). 

4. We neglect the effects of Doppler phase modulation due 
to LISA'S orbital motion which is an excellent approxi- 
mation for SMBH inspirals (S. A. Hu ghes, private com- 
munication, 2006: lArun et al.ll2007bl) . 

5. We neglect SMBH spins and, in particular, neglect the 
spin-orbit precession for determinations of the angular 
parameters. This assumption is useful in simplifying 
our equations and in focusing on the behavior of pure 
angular modulation. Including the effect of spin helps 
reduce uncer tainties in the final da ys and hours of the 
coalescence dLang & Hughes Il2007l) . 

6. We assume circular orbits (e = 0). Although eccen- 
tricity is effi ciently damp ed by gravitational radia- 
tion reaction dPetersI [19641) . the presence of gaseous 
circumbinary disks could lead to non-zero eccen- 
tricities fo r at least some LISA inspiral events 
jArmitage & Nataraian 2005; Pa paloizou et alj 120011 ; 
iMacFady en & Milosavljevic 2006j). For highly eccen- 
tric orbits, higher order harmonics appear in the GW 
phase, which could again lead to improved errors. 

7. We fix the start of LISA observations at a look-back 
time ti = min{f(/n,in), lyr} prior to merger This cor- 
responds to the time when the GW inspiral frequency / 
crosses the low frequency noise wall of the detector at 
/min = 0.03 mHz, but we limit the initial look-back time 
to a maximum of 1 yr before merger 

8. We assume that the slow annual amplitude modula- 
tion due to LISA'S orbital motion can be used to deter- 
mine the luminosity distance and angular parameters, 
Psiow = {dL,dN,(pN,dNL,<t'NL}, while the other remain- 
ing parameters, pfast = {A4:,M::fisco,</'isco}, are de- 
termined independently from the high frequency GW 
phase, with negligible cross-correlation between the 
two sets of parameters. This has been verified to be 
a good assumption up to a few days prior to ISCO, 

' The look-back time throughout this paper is defined to run "backwards" 
from ISCO, i.e. we set (merger = (see Paper I). 



with cross-correlations increasing the localization er- 
rors only by a few tens of %. Closer to merger, however, 
cross-correlations c an increase final e rrors at ISCO by 
a factor of 2-3 (La ng & Hughes Il2007h . 

9. We follow iBarack & Cutled (|2004 in calculating the 
spectral noise density, S„{f), which includes the in- 
strumental noise as well as galactic/extra-galactic back- 
grounds. For the instrumental noise, we use the 
effective non-angular averag ed onUne LIS A Sensitiv- 
ity Curve Generator"* (Bert i et alJ l2005ah . while we 
use the isotropic for mulae for the galactic a nd extra- 
galactic background dBarack & Cutled |2004[) . We do 
not c onsider errors from possib le theoretical uncertain- 
ties ( Vallisneri & Cutler1l2007h . which are negligible at 
look-back times ff well before fisco^ when the binary 
separations r » nsco- 

10. We estimate the expected parameter uncertainties us- 
ing the Fisher matrix approach, using the Harmonic 
Mode Decomposition (HMD) method developed in Pa- 
per I (essentially a discrete Fourier transform of the 
GW waveform with fundamental frequency yr"' cor- 
responding to the orbital period of the detector). The 
main advantage of the HMD method is that the time- 
dependence of errors can be extracted independently 
of the specific SMBH binary orbital elements, reduc- 
ing the computation time by many orders of magni- 
tude. This allows us to survey simultaneously the 
dependencies on source sky position, SMBH masses 
and redshift. We carry out Monte-Carlo (MC) cal- 
culations with 2x10"* random samples for the angles 
cos6n,cos9nl,(I)nl, and aj^isco- Here a is a time- 
independent combination of $(f) and S(f), defined as 
a = E-^ + (j)N-^. The angles are defined relative to 
LISA'S orbital plane (i.e., relative to the ecliptic), with 
6n and O^i referring to polar and (p^ and (j>NL referring 
to azimuthal coordinates. Several thousand values of M 
and z are considered, in the range 10^ < M/Mq < 10^ 
and 0.1 < z < 7. In addition, we ran specific MC calcu- 
lations to study possible systematic effects with respect 
to the source sky position (on a grid of several hundred 
values). 

2.2. Evolution of Error Ellipse 
2.2. 1 . Size of the Error Ellipse 

With the above assumptions, we obtain the two- 
dimensional sky position error ellipses. The time-evolving Icr 
errors on the size of the ellipse are shown for a few examples 
in Figure [T]'' We also note that the Fisher matrix approach, 
by construction, approximates the sky position errors (and all 
other parameter uncertainties) by ellipses. This approxima- 
tion has been shown for the 1 a uncertainties to be in excellent 
agreement with the posterior probability distributions for the 
slow parameters determin ed from Monte Carlo Markov Chain 
and B ayesian methods (iCornish & Ported 120061: IVallisneril 
l2007h . Figure [U shows the gradual improvement in locaUza- 
tion, for SMBH coalescences with six different fixed mass 
and redshift combinations, in the form of the major axis (2a), 
minor axis (2b) and equivalent diameter (2r = ^/4ab). 

^ www .srl.caltech.edu/~shane/sensitivity/ 
^ Note that throughout this paper, we show Icr errors. Higher confidence 
levels will obviously require covering proportionately larger areas. 
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Fig. 1 . — Evolution with pre-ISCO look-back time, tf, of LISA source localization errors for the sky position (major axis 2a, minor axis 2h, and equivalent 
diameter, 2r = 2\fah, of the error ellipsoid). Six cases with M = 2 X IO^-^ 'M© and z = 1,2 ai'e shown as labeled. Best, typical, and worst cases for random 
orientation events represent the 10%, 50%, and 90% levels of cumulative error distributions, respectively. The left panels show that the eccentricity of the angular 
error ellipsoid can change significantly during the last day. The absolute errors are typically small enough to fit in the FOV of a wide-field instrument in the last 
1-2 weeks. The horizontal lines show a diameter of 3.57°, which corresponds to localizing the source to within 10 deg^. 



Figure[T]displays results for three separate cumulative prob- 
ability distribution levels, 90%, 50%, 10%, so that 10% refers 
to the best 10% of all events, as sampled by the random dis- 
tribution of the five angular parameters. Note that the levels 
refer to the parameter in question (i.e. the 10% of all events 
for which a is best-measured, in general, is different from the 
10% of the events for which h is best measured), so eccentrici- 
ties can not be directly read off from this figure. The evolution 
of errors scales steeply with look-back time for tf > 40 days. 
For smaller look-back times, errors essentially stop improv- 
ing in the "worst" (90% level) case, improve with a relatively 
shallow slope for the "typical" (50% level) case, and improve 
more steeply in the "best" case (10% level among the real- 
izations of fiducial angular parameters). These behaviors are 
explained in Paper I; the minor-axis best case approaches the 
simple (S/N)~^ scaling expected in the absence of parameter 
degeneracies. 

In Figures |2] and |3] we show the probability distribution for 
the errors 2a and 2b in logarithmic bins at 4 different snap- 
shots: at ISCO, and at 1, 7, and 28 days before merger, us- 
ing 10^ random binaries. These distributions are broad, span- 
ning 1-2 orders of magnitude. They are generally asymmet- 
ric, skewed right or left for the major and minor axis, respec- 
tively. The distribution evolves with look-back time: it is ini- 
tially peaked, gets flatter and broader approaching merger, and 
develops a bimodal structure with two peaks. This bimodal 
nature of t he distribution at ISCO is also visible in the his- 
tograms of Lang & Hughes"! (12007 ). which our Figures |2] and 
[3] confirm at high resolution. The wide distributions suggest 
a non-negligible chance of being lucky /unlucky and finding a 
very well/poorly localized source. 



2.2.2. Eccentricity of the Error Ellipse 

Figures[ri|2]and|3]suggest that the 2D sky localization error 
ellipsoid is initially circular (a w b), but the geometry changes 
significantly during the last stages of the merger. For example, 
in the typical case (for z = 1 and mi = m2 =10*' Mq), in the last 
two weeks, the major axis improves much more slowly than 
the minor axis. The effect is much more pronounced at lower 
masses. As the left panels show, the major axis essentially 
stops improving within the last ^^day, while the minor axis 
maintains a steep evolution, so that the eccentricity of the 2D 
angular error ellipsoid changes rapidly (in fact this behavior is 
even more pronounced if spin precession effects are included, 
Lang & Hughes 2007). The fact that the error ellipse, in gen- 
eral, changes its eccentricity and its orientation may need to 
be taken into account in observational strategies (although, of 
course, as time evolves, subsequent error ellipsoids always lie 
entirely within the earlier error ellipsoids). 

Figure |4] shows the time-evolution of the axis ratio a/b 
of the 2D angular error ellipsoids for the same mass/redshift 
combinations as in Figures [T]-[3] To look for possible correla- 
tions between size and eccentricity, we first sorted all binaries 
according to the total area of the error ellipse into three bins, 
containing the "small", "medium", and "large" ellipsoids, cor- 
responding to 0-20%, 40-60%, and 80-100% of all events 
ranked by the instantaneous value of r = ^/ab. Figured then 
shows the evolution of the 10, 50, and 90 percentiles in ec- 
centricity in each of these three bins separately. The shaded 
region highlights the distribution of medium-sized ellipsoids, 
and the thick curves correspond to the median (50%) eccen- 
tricities in each bin. Depending on the binary parameters, we 
find a systematic trend between the size and eccentricity of 
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Fig. 2. — Probability density of the Ict LISA source localization ellipsoid major axis at 4 snapshots: ISCO, and 1, 7, and 28 days before merger Each curve 
shows the distribution in logarithmic bins, ln(2a), calculated using 10* random binaries. The vertical line shows the FOV of LSST for comparison. 
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Fig. 3. — Same as Figure|2]but for the sky locaUzation ellipsoid minor axis. 



the angular error ellipsoid. Figure|4]shows that if the ellipsoid 
is relatively small 3 months before merger, then it is very cir- 
cular, while the larger ellipsoids are more eccentric, but sub- 
sequently, at a few days before ISCO, the smaller ellipsoids 
become just as eccentric as the big ones. 

Overall, Figure |4] shows that the eccentricity evolution of 
the angular error ellipsoid is complicated, but that at the most 
relevant epochs, i.e. in the last two weeks prior to merger, 
when the absolute errors have shrunk to interesting levels, the 
eccentricity increases but remains small until the final day. 
More specifically. Figure |4] shows that the eccentricity is sig- 
nificant in the beginning, in the first few weeks after the binary 
enters LISA'S frequency band. However, the error ellipse then 



soon "circularizes", and the eccentricity remains small until 
the last day before merger The typical axis ratio decreases 
below 1.3, 1.1, and 1.01 after two months of observation for 
the large, medium and small ellipses, respectively. Except for 
the large ellipses, the eccentricity tends to increase again in 
the last week before merger The axis ratio is typically under 
1.3 before the last day to merger, with a 10% probability for 
a/b = 3. 

2.2.3. Effects of Mass, Mass Ratio, Upper Harmonics, and Spins 

For high values, the signal frequency at ISCO is low, 
and the binary coalesces soon after it enters LISA'S frequency 
band. For example, the inspiral of a binary with mi = m2 = 
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Fig. 4. — Evolution withpre-ISCO look-back time, tf, of the sky position error ellipsoid axis ratio (a/h) for the same six mass/redshift combinations (and equal 
masses = 1) as in Figure m Each group of three curves corresponds to the 10%, 50%, and 90% levels of eccentricity distributions, based on three bins (smallest, 
typical and largest) defined by the instantaneous value of the equivalent diameter r at each tf. Except for small observation times, the eccentricity remains small 
until the final day. This does not favor a strategy based on "tiling" the field with multiple pointings during a long-term monitoring campaign. 



10^ M0 and z > 1 can only be observed for less than 10 days. 
In such cases, the sky localization eccentricity does not have 
time to circularize and the eccentricity remains high (typically 
a/b ^ 5-10) throughout the entire observation. We note how- 
ever, that higher order GW harmonics beyond the restricted 
post-Newtonian approximation become relevant in this case, 
since these harmonics are inside the LISA frequency band 
much sooner. This improves angular errors gre atly for the 
high mass binaries jArun et alj|2007bt iTrias & Si ntes 2007). 
Due to the increase in observation time, we expect higher GW 
harmonics to circularize the angular ellipsoid for the more 
massive binaries.* 

Although Figures [T]-0] correspond to equal mass binaries, 
we also ran simulations for other mass ratios, q = m\/m2 ~ 
0. 1, and found very similar results. A plausible explanation is 
that two competing effects approximately cancel each other: 
(1) for fixed total mass M but lower mass ratio -q, the instan- 
taneous signal amplitude is lower {h+ x (0 « rj) while (2) the 
frequency evolves at a slower rate (d//df cx -q and the fre- 
quency at ISCO is independent of rf) implying that the binary 
spends more time at larger frequencies where the detector is 
more sensitive. Relative to the equal mass case, these com- 
peting effects approximately cancel out for ^ = 0.1 for both 
the size and eccentricity of the angular error ellipse. We find 
differences only for the high mass systems where the total ob- 
servation time is much larger for q = 0.\, leading to a smaller 
angular error ellipsoid and lower eccentricity. 

We expect this description to hold when including spin ef- 
fects until the final day prior to merger if the observation time 

' Note that the higher GW harmonics do not change the general expec- 
tation that the localization errors can be deduced from the slow amplitude 
modulation due to the detector motion. Therefore, compared to the restricted 
PN waveform, it is only the effective increase in the observation time which 
results in a difference for the localization ellipsoid. 



is at least months. For high spin binaries, however, spin pre- 
cession is expected to add an important extra feature to the 
waveform during the last day which improve s the final mea- 
surement angular erro rs by a factor of - 3 (IVecchio 1 12001 
iLang & Hug:hesll2006) . In this case, we expect spin preces- 
sion effects to further increase the angular e rror eccentricity 
during the last day. Indeed, recent results by Lang & Hughes I 
(2007) show that the minor axis acquires a long tail toward 
small values, 10"-1', indicating that the best events at ISCO 
will possibly be very elliptic. On the other hand, including 
spin precession reduces the error eccentricity compared to our 
estimates if the total observation time is less than a month (R. 
N. Lang, private communication, 2007). Therefore, spin pre- 
cession effects will have a tremendous impact on EM follow- 
up observations, allowing a 1-dimensional tiling of the field 
with small FOV instruments at late times. 

2.2.4. Sky position dependence 

Figures [T] and |4] show the basic statistics of the expected 
localization errors when all six angular parameters (2 sky po- 
sition, 2 angular momentum orientation, and 2 detector orien- 
tation parameters) are chosen randomly.^ This is useful to get 
an idea of the range of possible accuracies we may have when 
the LISA data becomes available. 

It is also interesting to calculate how the LISA localization 
errors depend on the actual sky position of the source, since 
EM counterpart detections can have very different prospects, 
depending on this position. For instance, it might be dif- 
ficult to find a counterpart if it lies in the Galactic plane, 
due to astrophysical foregrounds. On the other hand, a spe- 
cific example of a preferred direction could be the Lockman 

^ Note that the measured signal depends on only 5 comb inations of these 
angles, therefore only 5 variables are actually varied (§ 12. U . 
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hole, where the obscuring neutral hydrogen column density is 
anomalously low. 

To examine possible systematic effects with respect to 
the fiducial sky position, we have divided the (9^,(1)^) sky 
into cells and calculated the distribution of errors for each 
cell by choosing the remaining an gular parameters randomly . 
Surprisingly, in a previous work. lMoore & Hellingsl ( l2002h 
showed that no significant trends exist at ISCO as a func- 
tion of the polar angle 6*^. Our computations (for oti = 1112 = 
1O''M0 and lO^M© at z = 1) confirm this result to the 10% 
level, throughout the final 6 months of observation. (Median 
errors systematically get just slightly worse by less than 10% 
for sources that lie towards the Ecliptic.) Moreover, we find 
that localization uncertainties are independent of (p^. This, 
however, is not surprising since in Paper I we have shown that 
the LISA measurement depends on the combination - $, 
and is otherwise independent of ip^. Therefore the random- 
ness of $isco perfectly smear s out any systemat i c effec ts with 
respect to We note that iLang & Hughes I (l2007h find a 
very different 6^ dependence: they find Ecliptic equatorial 
sources to have ^ 2 times larger angular uncertainties than 
source s near the poles at ISCO. Although iLang & Hughes I 
(l2007h employ a different set of assumptions (accounting for 
correlations and spins, and choosing $ = S = at the start of 
the mission), these differences cannot be responsible for the 
strong ^^^^-dependence. The origin of this dependence there- 
fore remains unclear to us. 

On the other hand, it is plausible to expect systematic ef- 
fects in both the Ecliptic latitude and longitude if measuring 
the longitude relative to the detector position, (f)N-^{t). We 
find this to be a useful measure, since once the LISA plans 
are fixed, $(f) will be known, and 0^^- $(?) can be converted 
to Galactic coordinates for each (f)^ and t. For instance, de- 
pending on ^(January), it might turn out that a "good LISA 
direction" is in the Galactic plane in January and July, but not 
during other months. 

We calculate error distributions for cells of size AcosO^ x 
A((/)Ar - $isco) = 0.1 X 10° and find strong systematic ef- 
fects: the median localization error improves mostly in the 
instantaneous LISA plane during the final month, and it 
is worst in the two orthogonal directions ±A'^vo^st(Ol given 
by (6'iv,</>A.-$(f)) = (120°, 0°) and (60°, 180°). This re- 
sult has also been fo und independently in a recent study by 
iTrias & SintesI (l2007h using non-restricted PN waveforms for 
high mass binaries. The angular resolution is worst orthogo- 
nal to the LISA plane due to geometry: as our viewing angle 
of the binary changes, due to EISA's orbital and rolling mod- 
ulation, it causes the smallest change in the waveform since it 
corresponds to an extremum of the function [measured wave- 
form] vs [viewing angle]. Depending on the source position 
relative to the plane of the instantaneous orientation of the 
LISA triangle, the median angular localization uncertainty 
for r(f), in the case nii = 1112 = lO^M© and z = 1, ranges be- 
tween 5-50', 30'-60', 55'-90', 95'-135', and 170'-250' at 
ISCO, 1 day, 1 week, 1 month, 3 months, and 6 months be- 
fore ISCO, respectively. Therefore, during the most relevant 
last week, the area of the sky-position error ellipse can vary 
significantly, by a factor of 4-100, depending on source po- 
sition. We find that within the last month, the median errors 
are roughly rotation-symmetric around A'woist(f)-^ Therefore, 

* Note that this rotation symmetry might be a cons equen ce of choosing 
Sisco randomly in our Monte Carlo simulation (see § 12. It . We also tried 
fixing "3?isco — Sisco at some value and found an approximate rotation sym- 



systematic effects become prominent approaching ISCO be- 
cause errors improve much more quickly near the LISA plane. 
Quite remarkably, the errors are small enough very early on 
everywhere on the sky allowing to decide several months in 
advance whether the merger will take place in the sensitive 
region of LISA sky. Therefore, if several coincident inspiral 
events are identified during the LISA mission, it will be pos- 
sible to choose which one(s) to monitor with electromagnetic 
instruments according to which event is projected to have the 
best LISA precision closer to the merger. (The broader issue 
of cutting down the number of potential host galaxy counter- 
parts will be discussed in § l3.2.2l below.) 

3. TRIGGERED ELECTROMAGNETIC OBSERVATIONS 

The results summarized above suggest that it will be possi- 
ble to identify, prior to merger, a small enough region in the 
sky where any prompt electromagnetic (EM) counterpart to 
a LISA inspiral event will be located. Given sufficient "ad- 
vance notice," it will then be possible to trigger a search for 
EM counterparts as the merger proceeds (and also during the 
most energetic coalescence phase). 

The several square-degree field will, of course, contain a 
very large number of sources (a few x 10^ galaxies in total, at 
the limitin g optical magnitude of ^ 27 mag that may be rel- 
evant; e.g. iMadau | [T999 and discussion below). Having pre- 
dictions for the spectrum and time-dependence of a coalesc- 
ing SMBH binary, and therefore knowing what to look for, 
would obviously greatly help in identifying the counterpart, 
possibly allowing an identification even post-merger. Such 
a prediction, however, would require understanding the com- 
plex hydrodynamics and radiative properties of circumbinary 
gas at the relevant small separations of a few to a few thou- 
sand Schwarzschild radii. This is a notoriously difficult prob- 
lem even in the m uch simpler case of ste ady accretion onto a 
single SMBH (e.g.' Melia & Falcke II200T ). This suggests that 
the best strategy may be an "open-minded" search for any 
variable signatures prior and during coalescence. 

In the first half of this section, we describe several scenarios 
in which one may expect an EM signal to be associated wit h 
the coalescence of a SMBH binary (see also'Dotti et al.l2006l) . 
Most of this material in § 13. H is a review of existing literature. 
However, we also present some original material, such as cal- 
culations of the relevant orbital time-scales (Fig.|5]l, a rough 
estimate for the spin magnitude and orientation uncertainties 
for the binary prior to coalescence (Eq. [U, and we propose 
that candidates may be selected in advance based on whether 
or not the binary is expected to experience a large and appro- 
priately directed recoil that would favor bright EM emission 
(§13X2). 

In most cases, reliable assessments of the nature of the po- 
tential EM counterpart are difficult to obtain. What is most 
important for our purposes is that these scenarios make plau- 
sibility arguments for coalescing binaries being able to pro- 
duce detectable variable signals during their last few weeks 
prior to merger - either with periodic variations during the 
inspiral phase, or a strong transient signal during the final 
plunge/ri ng-d own/recoil phases. In the second half of this 
section, § 13.21 we give quantitative estimates of the capabili- 
ties of planned wide-field instruments to look for such vari- 
able signals. 

3.1. Expected or Plausible Variable EM Signatures 
metry in this case as well. 
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3.1.1. Periodic Flux Variations During the Inspiral Phase 

The problem of a SMBH binary embedded in an accre- 
tion disk has received relatively little attention in the lit- 
erature. It shares similarities with that of a planet em- 
bedded in the protoplanetary disk of its central star, which 
has been more extensively studied. The case of two near- 
equal mass SMBH, embedded in a thin, co-rotating disk, 
has been studied by a few authors. Numerical simulations 
jArtvmowicz & Lubow 1994VGunt her. Schafer & Klevl2004l: 
lEsc ala et al. 2004, 2005) indicate that, in this case, a central 
hole is expected to open in the disk due to tidal effects and 
the gas surface density within the binary's orbit could be sig- 
nificantly reduced. In the limit of efficient gap clearing, no 
EM signal would be expected until after a viscous time (t > 
years) past merger, when the inner edge of the gas disk falls 
in and starts accretin g onto the black hole remnant, p roducing 
an X-ray afterglow ( Milosavlievic & Phinnev Il2005h . 

However, in a more realistic situation, with no perfect az- 
imuthal symmetry, and especially for unequal masses or thick 
disks, some gas may still flow across the binary's orbit, and 
ultimately accrete onto one or both SMBHs. It is not un- 
reasonable to expect that this inflowing gas could produce 
a bright and variable source, with variability correlated with 
the orbital period. Indeed, the Newtonian potential around 
the binary will fluctuate periodically, with order unity varia- 
tions within a radius equal to the binary separation. The gas 
responds to this periodic perturbation, and a corresponding 
quasi-periodic flux variability h as been proposed on the basis 
of recent numerical simulations ( MacFadven & Milosavljevi3 
|2006; Hayasaki, Mineshige & Sudou 2007). 

Earlier numerical simulations, motivated largely to ex- 
plain the I2-yr variability of the blazar OJ 287, also found 
signifi cant mass in-flux across t he edge of the central 
hole dArtvmowicz & Lubowlll996l: iGiinther. Schafer & Kiev! 
|2004|) . especially when the disk gets thick (although these 
models relate to ~ 1000 times larger separations than the 
GW-emitting stage of interest here). Likewise, in the plan- 
etary context, close pre-main-sequence binary stars are ex- 
pected to cl ear central holes in their proto-planetary disks, but 
iJensen et~a i. ( 2007) argue that material can flow from the cir- 
cumbinary disk across the gap to explain the observed 19-day 
variability of the young binary UZ Tau E. 

There are other scenarios that may plausibly produce a 
variable signature related to the binary orbit. For example, 
lArmitage & Nataraian (2002) examine the case of a highly 
unequal mass binary, and argue that a significant amount of 
gas may be driven in by the merger if it is initially present 
interior to the orbit. Another way to promote variable sig- 
natures associated with the influence of the binary's orbital 
motion is to allow gas to flow in even as gravitational radia- 
tion rapidly shrinks the binary's orbit. This may be possible 
if the circumbinary disk is geometrically thick (or advection- 
dominated), with a viscous timescale much shorter than tradi- 
tional thin accretion disks. 

Assuming that variability is indeed present, what is the ex- 
pected period of this variability? Again, the answer is uncer- 
tain, but a reasonable guess is that the variability will be re- 
lated to half the orbital time-scale of the binary, correspond- 
ing to the periodic quadrupolar perturbation of the gravita- 
tional potential. 

In Figure |5] we show the time (At) spent by a binary with 
a total mass of lO^'^'^ M0, redshift z = 1, at each orbital pe- 
riod few (in logarithmic bins of few)- The plots show cases 
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Fig. 5. — The time spent (in hr) per logarithmic GW cycle-time bin during 
the inspiral for M = 2 X (10^ 10*, 1O')M0 (as labeled), q = 1 (sohd) and 
0.1 (dashed) for = 1 in the restricted 2.5PN approximation. In each case 
two curves with the same style correspond to the two extremes: maximum 
aligned and anti-aligned spins; a general case lies in between these extremes. 
Tgw is the GW cycle time in hours. The curves end at fj = 1 yr or at /„,i„ = 3 X 
10"' Hz, coiTesponding to Tqw = 9.3 hr. Note that the number of GW cycles 
can be estimated as A'cycics = y/-^, from the x and y-axis values corresponding 
to each curve. 

with mass ratios of ^ = 1 (solid) and 0.1 (dashed). The vari- 
ous curves of the same type show the upper and lower bound 
of the time-evolution depending on the magnitude and orien- 
tation of BH spins. The presence of spin modifies appreciably 
only the final day of the merger 

Quite remarkably. Figure |5] shows that the expected EM 
variability timescale is in a very fortunate range for detection: 
between several minutes up to 10 hours within the last year 
of merger. Detecting slower variations would require several 
exposures over exceedingly long time-periods, while detect- 
ing faster variations (of the same amplitude) would require 
much higher sensitivity per unit exposure time. If a periodic 
electromagnetic variation is found within this range, e.g. for 
an AGN in X-rays, it would indicate evidence for an SMBH 
binary approaching coalescence. The evidence would be es- 
pecially convincing if this nearly periodic variable "light- 
curve" tracked the known, nearly periodic GW emission, with 
a significant cross-correlation between the periods of the two 
waveforms. The EM lightcurve may perhaps even follow the 
slow decrease in period, related to the GW chirp, shown in 
Figure |5] 

On the basis of the relatively well-defined nature of these 
potential counterparts, it may even be interesting to carry out, 
well before the LISA mission becomes operational, system- 
atic searches for such quasi-periodic behavior with the various 
all sky variability surveys currently being planned. 

3.1.2. Transient Signatures During Final Coalescence 

Perhaps the most promising phase for the emergence of 
bright EM counterparts to SMBH merger events is, in fact, 
during or shortly after the violent final coalescence stage. 
This is the most energetic phase of the binary's evolution, 
with typically ~ 10^** ergs of gravitational radiation released 
to the environment in a short ^ 10^ s. There are at least 
two ways to couple this energetic outburst to the imme- 
diate gaseous environment o f the SMBH binary (see also 
iMilosavIievic & Phinn ev 2005). 

A few percent of the binary system's rest mass is lost when 
the binding energy of the two holes is carried away by GWs. 



9 



Most of this loss occurs during the final cycles of inspiral, 
the plunge and subsequent ring-down. Gas bound to the 
SMBH's binary would thus experience a sudden reduction in 
the mass of the central object, on a timescale typically much 
shorter than the longer (local) dynamical timescale on which 
this gas can respond. Although this sudden disturbance is 
of moderate amplitude a few %), it is a potential source 
of transient variability associated with the coalescence, de- 
pending on details of th e disturbed surrounding gas properties 
( iBode & Phinne v 2007). 

A second, and possibly more direct and promising avenue 
to couple the energetic GW burst at coalescence to the SMBH 
binary's environment is through the GW recoil/kick experi- 
enced by the system. The accumulation of linear momentum 
leading to the final recoil builds during the final cycles of in- 
spiral, the plunge and subsequent ring-down, so that it appears 
as quasi-instantaneous except for any gas that may be located 
in the immediate vicinity of the coalescing binary. The magni- 
tude of the kick is determined by the mass ratio, the magnitude 
of the spins, and the relative orientation of the spin and orbital 
angular momentum vectors. The recoil kicks can be calcu- 
lated accurately using numerical relativity (e.g. Ba ker et alj 
I2OO6I) . The final BH velocity is found to reach 170km/s 
for no spins dGonzalez et alj|2007h . a nd 2500-4000 km/s for 
high spins in special configurations ('Campanelli et al.ll2007l; 
iBriigmann et al. 2007). For comparable mass binaries with 
large spins and arbit rary orientation, a recoil > 500 k m/s has 
a 30% probability dSchnittman & Buonanno1l2007h . If the 
spin axis is aligned with the angular momentum, as maybe 
the case in gas rich mergers, then the recoil is expected to be 
< 200 km/s (Bogdanovic et al. 2007, and see § l3.2.2l below). 
Even for a modest recoil velocity of 100 km/ s, this amounts to 
^ 10''^ ergs of mechanical energy that can in principle be de- 
posited in the environment as a SMBH remnant of ^ 10^ M0 
ploughs through. The EM signatures that may emerge from 
such interactions are difficult to evaluate a priori since they 
will depend on many details, such as the properties of the gas 
that remains bound to the recoiling SMBH, the geometry of 
the recoil system with respect to its environment, and possibly 
the remnant BH spin. We will simply note here that these EM 
counterparts are not necessarily bound by the Eddington limit 
traditionally applied to steady systems, if the recoil energy is 
deposited fast and efficiently enough for the event to adopt an 
outflowing/explosive character. 

For the present discussion, it is interesting to ask whether 
the final recoil magnitude and direction can be estimated in 
advance, and to determine whether it lies in the plane of the 
binary, where a gaseous disk may be preferentially located. 
If so, then this advance knowledge could possibly be used 
to pick the events that are most likely to produce a transient 
EM counterpart. Since recoil velocities are determined by the 
mass ratio, spin magnitudes and orientations, the question is 
whether these parameters can be estimated with sufficient pre- 
cision in advance. Without specific computations, we can es- 
timate the advance uncertainties for these parameters using 
simple scaling arguments. The mass ratio is determined from 
the IPN effect, hence it is a "fast parameter" that will likely 
improve quickly with time, close to (S/N)~^ (see § 12. 11 1. Spin 
effects arise at 1.5PN order, implying that the binary spins 
are also "fast parameters". The final errors on the spin mag- 
nitudes are excellent (see ,Lang & Hughes 2006), better than 
l%-6%-20% for z = 1-3-5 if mi = ma = IO^'Mq. Even the 
final spin angle uncertainties have not been estimated previ- 



ously, but we can assume that they are not much worse for 
large spins than the angular momentum orientation uncertain- 
ties {S9l, S(j)L), since the two are coupled through precession, 
and precession was shown to help in reducing the {6di,S4>L) 
errors. If the spin magnitude is chosen randomly, then the spin 
orientation measurement error is 3.5 deg for mi = m2 = lO^M© 
at z = 1 (R. N. Lang, private communication, 2007). We as- 
sume conservatively that the evolution of spin errors scale 
with (S/N)~^ until the last day before ISCO. According to 

1 /3 

Paper I, S /N oc t^' .In this case 
/ tf 

5spm(f/)~ 5spm(fisco), (1) 

V 'ISCO / 

with risco = 33minx(Mj/(4 x lO^M©)). For mi = m2 = 
IO^'Mq at z = 1, we find ^spin ~ 1% x [f//(33min)]i/l i.e. 
spin magnitudes will be determined to ~ 8% at 10 day before 
ISCO. For the direction, ^10 deg 20 deg) precision can be 
achieved 10 days before ISCO, if the final uncertainty at ISCO 
is 1.3 deg (2.6 deg). This information could thus in principle 
be used to focus on events with specific magnitudes and/or di- 
rections for the expected recoil, in an effort to maximize the 
likelihood of finding bright EM counterparts. 

3.2. Searching for Variable Signals with Wide Field 
Instruments 

3.2.1. Basic Parameters for a Variability Search 

The most obvious observational strategy to catch a prompt 
EM counterpart, especially in the absence of robust theoret- 
ical predictions, is to blindly and continuously monitor the 
source area, through coalescence, for variable sources. The 
requirements for such a monitoring - large FOV, and a fast 
camera - coincide with those of searches for distant Super- 
novae (SNe) in the optical bands. Motivated by the latter, 
there are several instruments being built or planned with spec- 
ifications that could prove suitable for a LISA counterpart 
search. Table[T]shows a non-exhaustive list of facilities, either 
planned or being built, that have a several square-degree field 
of view and may be suitable for a LISA counterpart search. 

According to the arguments in the previous section, a co- 
alescing binary may produce some variable signature in the 
last two weeks before ISCO. The amplitude and period of any 
variability, however is highly uncertain. In view of this un- 
certainty, it is interesting to work backwards, and compute 
the expected brightness of a variable signal that is some fixed 
fraction of the Eddington luminosity. An Eddington-luminous 
SMBH at z = 1-3, with a mass of (1O''-1O^)M0, assuming a 
10% bolometric correction, would appear as a ~ 22-25 mag 
point source in the optical bands. This implies, for exam- 
ple, that if the variable signal is ^ 10% of the Eddington 
luminosity, then we need sensitivities corresponding to 24- 
27 mag. Ten times lower mass BHs would require sensitivities 
of 26-29 mag. This variability is to be detected on an orbital 
time-scale, which is between several minutes to several hours 
(Fig.H). 

Scaling from the sensitivities shown in Table [T] assuming 
S/N (X \fi, we see that a sensitivity of 25-27 mag can be 
reached in ^ Imin to Ihr integrations with Pan-STARRS 
and with LSST We conclude that an LSST-class instrument, 
repeatedly integrating over the LISA area, would be able to de- 
tect a variability that corresponds to ^ 10% of the Eddington 
luminosity and follows the period of the binary's orbit. 

Figures [T] and |4] show that the eccentricity of the LISA error 
box may be significant during the last day or for high-mass 
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TABLE 1 

Wide Field Electromagnetic Instruments 



Name 


FOV (deg-) 


bands 


mag. limit (exp) 


start date 


reference 


High Energy Transient Explorer-2 


0.9, 1.6 srad 


soft, hard X-ray, 7 


8x10 — S- 

cm^s 


2000 


space.mit.edu/HETE 


Galaxy Evolution Explorer 


1.2 


Near, Far UV 


21-22 (1 hr) 


2003 


www.galex.caltech.edu 


Palomar-Quest 


3.6 X 4.6 


U,B,R,I,r,i,z 


17.5-21.0 (140sec) 


2003 


B_altav et al. (2007) 


XBootes Survey 


9.3 


X-ray 


(l-8)x 10-'5^ 


2005 


www.noao.edu/noao/noaodeep 












/XBootesPublic 


SkyMapper 


2.3 X 2.4 


u,v,g,r,i,z. 


20.6-21.9 (llOsec) 


2007 


w w w. mso . anu . edu . au/sky mapper 


South Pole Submillimeter Telescope 


1 


sub-mm (350-850)/im 


1 mJy (18-27 hr) 


2007 


cfa-www.harvard.edu/~aas/tenmeter 


Low Frequency Array 


20 


radio (10-240)MHz 


(0.03-2)mJy (1 hr) 


2007 


www.lofar.org 


Pan-STARRS-1 


2.5 X 2.5 


g,r,i,z 


24 (1 min, 5a-) 


2008 


pan-starrs.ifa.hawaii.edu 


Gamma-ray Large Area Space Tel. 


2 srad 


7 


6x 10"'-!^ (lyr) 


2008 


glast.gsfc.nasa.gov 


Dark Energy Survey 


2.2x2.2 


g,r,i,z 


24(100min,5cr) 


2009 


www.darkenergysurvey.org 


Kepler 


105 


430-890nm 


24(150min,5cr) 


2009 


kepler.nasa.gov 


Large Synoptic Survey Telescope 


10 


u,g,rJ,z,Y 


27(lhr, 5cr) 


2014 


www.lsst.org 


Joint Dark Energy Mission 


15 


6 optical -1- 3 NIR 


27 (1 hr) 


2016 


universe.nasa.gov/program 












/probes/jdem.html 


Wide Field Imager 


1 


350-1800 nm 


24 (110 min) 


2017 


sci.esa.int/science-e/www/object 








2x 10-"^ 

cm'-s 




/index. cfm?fobjectid=39692 


Energetic X-ray Imaging Survey Tel. 


154x65 


X-ray 


2020 


exist.gsfc.nasa.gov 


Square Kilometer Array 


1 


radio (20-200)GHz 


1 AtJy (1 hr) 


2020 


www.skatelescope.org 



Note. — The table shows EM instruments, existing, being built, or planned, that have a several square-degree field of view and may be suitable for 
a LISA counterpart search. The 3rd column shows the observational wavelength bands, and the 4th column shows the sensitivity that can be reached in 
these bands in the indicated integration time at the quoted S/N (derived from information at the websites listed for each instrument in the 6th column, when 
possible). 



SMBHs with M > 2 x IO^Mq, with the more uncertain (ma- 
jor) axis longer by up to a factor of several than the better 
determined (minor) axis. This may impact the plausible du- 
ration of any triggered monitoring. If the telescope is rapid 
enough, then one requires only that the minor axis fits within 
the FOV (assuming that the FOV is spherical), since it will 
be possible to tile the rest of the LISA ellipse by a few ex- 
tra pointed observations, displaced along the major axis. On 
the other hand, if the required integration time is comparable 
to the expected orbital period, then performing such a tiling 
would mean shallow observations that may miss the variabil- 
ity. In this case, the monitoring campaign can start only after 
the major axis has shrunk below the linear size of the FOV. 

Figures |6]and|7]display advance warning time contours for 
typical (50%) and best case (10%) events, adopting the LSST 
FOV as a reference. Advance warning time contours are log- 
arithmically spaced, with solid contours every decade and the 
shaded region highlights the (M,z) region where at least 10- 
day advance notice will be available. The top panel of Fig- 
ure |6] shows that 10 day advance warning to cover the full er- 
ror ellipsoid with a single LSST pointing is possible for a large 
range of masses and source redshifts, up to M ^ 3 x lO^M© 
and 1 .7. The bottom panel shows that the (M,z) region in 
which this same 10 day warning is limited to the minor axis 
of the error ellipsoid reaches to somewhat higher redshifts. 
Figure |7] shows how far the advance warning concept can be 
stretched, by focusing on the 10% best cases of random ori- 
entation events. The top panel shows that a 10 day warning to 
cover the full error ellipsoid with a single LSST-type pointing 
is possible up to z 3 for masses around M ^ 1O*M0. The 
bottom panel shows that the (M,z) region in which the same 
10 day warning is restricted to just the minor axis of the er- 
ror ellipsoid reaches out to comparable redshifts, for a similar 
range of masses. 

As explained above, a unique variable EM signal could 
be produced in the last stages of the merger. Figures |6] 
and Figure [T] show that requiring a warning of just one day 
would extend considerably the range of masses and redshifts 
for which a single LSST-type pointing is sufficient, out to 



z > 5 for the bes t events. However, as we mentioned above, 
iLang & Hughes I (l2007h find that in the last day, parameter 
correlations can degrade the errors by a factor of 2-3, but this 
degradation is more than offset by the factor of several im- 
provement in the errors caused by spin precession (for high 
spin). Therefore, our localization errors for the last day are 
overly optimistic, by a factor of 2-3, for low spin mergers, 
but are conservative for high-spin mergers. Also, the issue of 
eccentricity turns out to be most important in this regime (i.e. 
for localizing the be st even ts ^ 1 day prior to ISCO, especially 
for high spins, see § l2.2.3l l. In this case, requiring the FOV to 
cover only the minor axis would allow going out to z ~ 8. 

It should be possible to significantly improve on the above 
magnitude requirements, if the signal is indeed close to peri- 
odic at few- This is because the above calculation requires the 
flux variations to be detectable for any individual source in the 
LSST FOV. A lower-amplitude variable signal would be lost 
in the noise for any individual source (i.e. it will be consistent 
with photometric noise). Fitting to a sinusoid-like template 
would require much less sensitivity. Even in the absence of 
a precise template for the EM emission, one could use the 
GW waveform itself (or variations of it). Cross-correlating 
the EM flux with the peaks and troughs of the GW waveform 
should allow one to recover such below-the-threshold vari- 
ability, effectively decreasing th e small est detectable variabil- 
ity amplitude by a factor of ^ ^A'cycies, where A'cycies ~ lO-''^'' 
is the number of orbital cycles for the binary during the last 
two weeks, forM = 1O'''*^M0, respectively (see Fig.|5]l.^ This 
type of template matching could presumably be attempted in 
various regions of the EM spectrum (see Table 1 for a few 
possibilities). 

3.2.2. Cutting Down on the Counterpart Candidates 

The above scenario so far envisions monitoring all sources 
in a few square degrees, in ^ hourly intervals, for ^ two 
weeks prior to and around ISCO. This wide field will contain a 

' Note that ^A'cycies scaling is connect only for constant modulation power 
per cycle, which might not be the case for the EM variations. 
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Fig. 6. — Contours of advance warning times in the total mass (M) and 
redshift (z) plane, for typical events (50% level of cumulative distributions for 
random orientation events) with SMBH mass ratio mi/m^ = 1. The contours 
trace the look-back times at which the major axis {2a, top panel) and the 
minor axis (2h, bottom panel) of the localization error ellipsoid first reach 
an LSST-equivalent field-of-view (3.57°). The contours are logarithmically 
spaced in days and 10 days is highlighted with a thick curve. 

very large number of galaxies. As seen earlier, the interesting 
optical magnitude limit may be ~ 25-27, for variability at the 
level of 1 - 10% of the Eddington limit. The surface density 
of galaxies at this magnitude is known from deep op tical ob- 
servations to be about ^ 10^galaxies/(10deg^) (e.g. lMadaul 
[T999.) . While this large number of possible candidates may 
seem daunting, there are several ways to significantly shrink 
the candidate list. We enumerate these options in this section. 

Final Angular Localization Cut. In Figures |6] and |7] we 
showed the advance warning times when the source localiza- 
tion can be first covered with a lOdeg^ instrument. However 
not all galaxies are going to be relevant in this region even- 
tually, since the area of the error ellipsoid shrinks continu- 
ously (Fig. [T]i. For m\ =m2 = IO^Mq at z = 1, for instance, 
the lOdeg^ localization is typically reached at ff w 86 days. 
The (10-90)% distribution levels of the sky position uncer- 
tainty later decrease to 5Vl = (0.1-1) and (0.6-7) deg^ at ff = 1 
and 10 days before ISCO (Fig. During the last day, the 
size of the ellipsoid decr eases by another factor of ^ 3 for 
a source with large spins (iLang & Hughes Il2007h . Therefore 
the number of candidates among the galaxies within the orig- 
inal lOdeg^ field decreases to a fraction (6-70)%, (1-10)%, 
and (0.3-10)%, at 10 days, and 1 day before merger, and at 
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Fig. 7. — Same as Figure|6] except that results for the best (10% cumulative 
level) events are shown. 

ISCO, respectively. Obviously, these cuts can only be applied 
progressively as the merger progresses, always abandoning 
candidates that move out of the shrinking error ellipse, and 
not immediately at the beginning of the monitoring campaign. 

Photometric Redshift Slices. All of the survey instruments 
listed in Table 1 should be able to observe the field in several 
bands, and, for example, LSST should be able deliver photo- 
metric redshifts of all sources with an accuracy of 6z < 0.05 
(see, e.g. iHuterer et alj[200 6'). If LSST has already surveyed 
the target area before the LISA event, then galaxies will have 
tabulated magnitudes and photo-z's already. These data may, 
however, not pre-exist, or may not go deep enough (if the 
SMBH system has too small a mass). But the cumulative ex- 
posures of this particular deg^ field during a monitoring 
campaign should be sufficient to reach target errors already in 
approximately the first day of monitoring. 

Importantly, LISA will deliver the luminosity distance diiz) 
with high accuracy. In Figure [8] we show the errors on dis- 
tance as a function of time. For a flat ACDM cosmology with 
(n\,ilm,h) = (0.7,0.3,0.65) the redshift uncertainty is given 
by & = (0.8 1 , 1 .6)6di ^/dL foi";^ = (1; 2), respectively, provided 
that Sz > 0.005 ( Kocsis et al.ll2006h . Few percent errors are 
thus reached several days prior to ISCO, except for the worst- 
case events, in the entire mass range, 10^-10^ M0, out to 
z = 2. However, the distance determination is ultimately lim- 
ited by weak lensing uncertainties due to fluctuations in the 
dark matter density along the line of sight, shown as a hor- 
izontal line in Fig g] at &wi = (0.025,0.097) for z = (1,2), 
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that will host the SMBHs with Msmbh = IO^'^Mr, a conser- 
vative assumption is S/T = 0.1-1 (Benson et al. 2007). Inte- 
grating the observ ed luminosity func tion of galaxies between 
these Mb bounds dRvan et al.ll2007l using their fits for faint 
galaxies at z = 1, $* = 26.1 x IQ-^U^c'^ , M| = -22 4, a(z) = 
-l.12 -0.12z, and P = -1.00, g = 1.72 given bv iLin et all 
|1999|) for ^0 = 0.1, and extrapolating this relationship to low 
luminositie s) and using th e cosmological comoving volume 
element dEisenstein Ill997h 



Fig. 8. — Evolution with pre-ISCO look-back time, tf, of LISA errors on the luminosity distance dtiz) for mi = m2 = (10^ IC*, 10') M© and z = (1,2) 
as labeled. The 10%, 50%, and 90% levels of cumulative error distributions for random orientation events are shown. Horizontal lines delineate the level of 
uncoiTected weak lensing uncertainties. 

respectively. The weak lensing uncertainty cannot e asily be 
corrected by more than a fact or of ~ 40% (see Da lai et alj 
|2003; Kocsis et al. 2006; Gun narsson et al.ll2006l) . Once the 
LISA uncertainty decreases below this limit, the total redshift 
error will not improve. The redshift distribution of optical 
galaxies at 27-2 8 mag is already known from deep observa- 
tions (e.g lCoe e t al. 2006, using the Hubble Ultra Deep Field), 
with roughly 1-10% of them falling in a radial shell of width 
Sz ~ 0.05 around redshift z = 1-2, hence photometric redshifts 
will offer a factor of ^ 10-100 reduction in the number of 
candidates. 

Expected Luminosity Range. Since LISA will also mea- 
sure accurate BH masse s prior to the merger (the masses are 
"fast parameters", see § 13.1.2b it may be possible to roughly 
estimate the magnitude of the host galaxy from the known 
empirical correlations between galaxy properties and SMBH 
mass. For example, extrapolating from the sample of nuclear 
SMBHs in the local universe, a lO*" M© SMBH is expected 
to reside in a galaxy with a spheroid mass ~ 10^ Mp^, o r abso- 
lute magnitude Mb ^ -16 (iFerrarese & Merritt1l2000h . corre- 
sponding to an apparent optical magnitude of m '-^ 28 at z = 1 . 

A large fraction of the galaxies will not have luminosities 
consistent with the mass-luminosity relationship and the BH 
masses determined by LISA, allowing to ef ficiently reduce the 
numbe r of counterparts. Using Eq. (19) of IFerrarese & Ford] 
(12005b . we find that the B-band spheroidal luminosities 
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10^ X (7.8,3.7, 1.7) and 
10^ X (10,3.6, 1.2) for Msmbh = 2 x IO'^^^-^'Mq, Sz = 0.05 and 
for z = 1 and 2, respectively. If a bulge/disk decomposition 
is also available to constrain S/T (Ben son et al.i 120071) . then 
A^coiinterpait is rcduccd by another factor of ^ 2. Therefore 
imposing bounds on the range of both the luminosity and the 
redshift of the counterpart reduces the number of counterparts 
significantly. Incorporating the angular localization cut above 



leads to a total of Nr 



counterpart 



areMB = (-12.7!^;^,-15.l!| 



1+1.6 

.1 ' 



17.5!^:°,) for SMBH masses 
A^SMBH = 2 X 10^'"' 'Mq. Note, that the eiTors are (by far) dom- 
inated by the uncertainties in the M^mbh-Mb relationship, and 
not the LISA measurement errors, with (5Msmbh/A^smbh < 
10"'*. In order to obtain the total luminosity, we must divide 
the spheroidal luminosity with the spheroidal/total luminosity 
ratio S/T. This ratio depends on the galaxy morphology, ori- 
entation, and luminosity. For the low-luminosity spheroids 



lO'-lO'* at 10 days before 
merger, 100-1000 at 1 day, and 1-1000 at ISCO. 

We caution, however, that before applying cuts on candi- 
dates, one should be well aware that the relation between to- 
tal SMBH mass and spheroidal luminosity is currently known 
only for nearby galaxies, selected by observational feasibility 
criteria. Clearly, the relation can be different for rare objects 
of a transitional nature, such as SMBHs in the last few weeks 
of their merger 

Expected velocity dispersion. There is observational evi- 
dence for a tight correlation between the SMBH mass, M, 
and the large scale bu lge velocity dispersion of galaxies, a 
dGebhardt et al.l 120001) . By contrasting M, as measured by 
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LISA, with a obtained from galactic spectra, a major frac- 
tion of the gal axy candidates could be ruled out. Using 
equation 20 in iFerrarese & Ford I (l2005h . and extrapolating 
to lower masses, we find a = (50 ± 5, 81 ± 4, 130 ± 1), for 
M = 2 X 10^ '''^Mq, respectively. In comparison, the veloc- 
ity dispersi on for observed ga laxies ranges between a ^ (60- 
400)km/s dSheth et al.ll2003h . Since we expect there will be 
many more galaxies with low velocity dispersion, a very effi- 
cient reduction of A'counterpan will be only possible for the most 
massive detected SMBHs. Unfortunately, the measurement 
of velocity dispersions for faint galaxies is extremely difficult 
in practice. It may be feasible only for the larger mass BHs, 
and presumably after longer-integration spectra are obtained 
well after the merger As before, we caution that the rela- 
tion between total SMBH mass and the velocity dispersion 
is unknown for rare objects of a transitional nature, such as 
SMBHs in the process of coalescence. 

Quasar counterparts. If the LISA source is producing a 
continuous near-Eddington luminosity with fluctuations at 
the few percent level, then most likely, it will appear as a 
point source (i.e. a faint quasar) and the host galaxy will not 
be detectable. The bright point source in the nucleus, with 
^ 24 mag, will outshine the galaxy. Kocsis et al. (2006) found 
that the typical LISA error box {5^ = 0.3 deg^, & = 2.5%), 
for a merger event produced by ^ 10^ SMBHs at z = 1 
would statistically contain ^ 1 quasar with at least 10% Ed- 
dington luminosity; the ^ 30 times larger area ~ 2 weeks be- 
fore merger will contain ~ 30 such quasars. These quasars 
could be identified in advance, and monitored separately for 
variability with smaller-field of view instruments. 

Orbital Eccentricity as Indication of Circumbin ary Gas. In 
2D sim ulat ions of circumbinary disks, ^ mitage & NataraianI 
(120021) and iMacFadven & Milosavhevid (i2006) find that the 
interaction of the binary with the ambient gas produces ec- 
centricity (in the orbit of both the gas and the binary). The 
eccentricity leaves a s t rong impri nt on the GW waveform 
(iBarack & Cutleiil2004l: IVasuth & M ajar 2007). Moreover, it 
is a "fast parameter" since it modifies the spectrum of the sig- 
nal on the binary orbital timescale (see § 12. lb . Therefore we 
expect that the LISA data-stream should be able to discover 
such eccentricity several months before merger and thus alert 
us to the occurrence of a gas-driven inspiral. 

Spin and Orbital Plane Alignments as Indication of Am- 
bient Gas. Prior to the merger, the two SMBHs may reside 
in two bulges, which are driven together during the galaxy 
merger by large-scale torques (e.g., Dubinski et al. 1996). As 
a result, the orbital plane of the binary SMBH may be aligned 
with the plane of the galactic remnant. Indeed, in gas rich 
mergers, the binary plane may align with the plane of the disk 
dBogdanovic et alj|2007h . Moreover, these authors argue that 
the interaction with the gas also aligns the individual BH spins 
in the direction of the orbital angular momentum. Therefore 
the relative orientation of the binary orbit, spin vectors, and 
embedding galaxy may be valuable observables. In particu- 
lar, if the measured spin vectors are found to align with the 
orbital angular momentum, it may indicate that a sufficient 
amount of gas was present in the history of the merger. More- 
over, in the case of alignment, the direction of the recoil kick 
at merger is in the orbital plane, making the kick-triggered 
variability discussed in § 13.1.21 more plausible. Conversely, 
misaligned spins may indicate that the merger is dry, and an 
EM signal might not as easily accompany the GW event. Note 
that a ±20 deg triple alignment between the spins and orbital 



momentum by chance has a mere 9 x 10""* probability. In 
case there are several coincident inspirals present in the LISA 
data-stream, a favorable EM follow-up candidate could thus 
be chosen according to a criterion for potential EM activity 
based on the observed relative alignments. 

In Figure|9] we show the accuracy with which we will know 
the inclination of the SMBH orbital angular momentum, in 
advance of the merger, as specified by the two angles {fii, (pi). 
When considering the individual EM candidates in the LISA 
error volume, one-by-one, the EM information allows to set 
the sky position errors {&9nt54)n) to essentially negligible val- 
ues. Therefore, it is meaningful to calculate the 2D error el- 
lipsoid for {9l,4>l) by neglecting correlations with ddj^ and 
6(pN- Figure |9] then shows that we will know the inclination 
remarkably well in advance. One may therefore envision fo- 
cusing on galaxy candidates that exhibit inclinations similar to 
that of the SMBH binary. This will be feasible only if the host 
galaxy is actually detected and resolved (i.e. if the SMBH bi- 
nary does not produce near-Eddington emission whose PSF 
completely hides the galaxy). The angular size of small el- 
liptical galaxies at z = 1-3 with optical magnitudes ~ 27 is 
expected to be 0. 1-0.5 arcseconds (see, e.g., simulations by 
iKawata & Gibsonll2003h . which cannot be resolved from the 
ground. However, by z ~ 1, most of the host spheroids may 
be surrounded by much larger disks. For example, the 20 kpc 
diameter disk of a galaxy like our own Milky Way would ex- 
tend beyond 2 arcseconds at z = 1-3, and could thus be re- 
solved. Again, one can envision separately monitoring with 
smaller field-of-view instruments a few of those galaxies (if 
any) whose disks appear to be closely aligned with the orbital 
plane of the SMBH binary. 

In conclusion, the alignment of spins (indicative of the pres- 
ence of gas) and the alignment of orbital momentum with the 
plane of the galaxy can be very helpful in choosing the most 
likely host galaxy candidates. Note, that a ±20deg coinci- 
dence between the GW and EM measured orbital orientation 
has a 3% probability, thereby allowing a very effective cut to 
be made on the candidates. 

Galaxy Overdensity as Indication of Enhanced Merger 
Rate. Another possibility is to focus on galaxy clusters 
within the few square degrees of interest. There may be 
a few to a few tens of known galaxy clusters in this area. 
Combining with the photometric redshift information, Az = 
0.05, there are on average A^ciuster/ deg^ = 0.2 and 0.01 clus- 
ters within this shell, a t z = 1 and 2, respectively (e.g. 
Holder, Haiman, & Mohr ^ 200 Ih . The probability that a ran- 
domly chosen galaxy has undergone a recent merger is sig- 
nificantly larger in a cluster than in the field; again, it would 
make sense to separately monitor the dense cluster environ- 
ments with smaller FOV telescopes. 

4. MOTIVATIONS AND NEW FUNDAMENTAL TESTS OF GRAVITY 

We have argued that real-time triggered searches for EM 
counterparts to SMBH binary merger events observed by 
LISA will be technically feasible, at least for a subset of all 
such events. Observational strategies to achieve such identifi- 
cations constitute major astronomical endeavors, however It 
is thus important to clarify the nature of the additional sci- 
ence that could be enabled by real-time identifications. We 
describe two important motivations for pre-merger localiza- 
tions here. 

The first motivation is simply coincidence. Post-merger 
identification of host galaxies will likely rely on active nu- 
cleus variability or on a statistical association with rare galac- 
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Fig. 9. — Evolution with pre-ISCO look-back time, tf, of LISA errors for the orientation of the orbital plane of the SMBH binary relative to the Hne-of-sight. 
The inclination and rotation uncertainties are correlated leading to an error ellipsoid with major axis, 2aNL, equivalent diameter 2rNL and minor axis Ib^L. Best, 
typical, and worst cases for random orientation events represent the 10%, 50%, and 90% levels of cumulative error distributions, respectively. The range spanned 
by the Ir^L distribution is highlighted for clarity. The figure shows that the binary inclination angle can be precisely localized very early on, to within degrees, 
allowing a selection of prime host galaxy candidates that have a similar inclination. 



tic objects (e.g. quasars, ULIRGs, X-ray or starburst galax- 
ies; iKocsis etal . 2006; Dotti et al. 2006). Depending on the 
nature of the AGN variability or the specific galactic attribute 
one is willing to associate with a SMBH binary merger, pos- 
sibilities for chance coincidences may remain uncomfortably 
high. The combination of space (sky location) and time coin- 
cidence for a prompt EM counterpart that is gravitationally- 
timed with a SMBH merger should considerably reduce the 
risks of chance coincidences and would thus likely be the 
most secure way of identifying a truly unique counterpart 
and host galaxy. This is particularly important in this con- 
text, with observations of a rare transient astrophysical phe- 
nomenon that has never been observed before, whether grav- 
itationally or electromagnetically. Secure counterpart identi- 
fications are also crucial if one is to perform new fundamen- 
tal tests of gravitational physics with these events, such as 
comparisons bet ween electromagnetic a nd gravitational Hub- 
ble diagrams (Deffayet & Men oullIOOTi) . 

The second important motivation is that the detection of 
prompt electromagnetic counterparts itself may provide ad- 
ditional opportunities for testing gravitational physics on cos- 
mological scales. We highlight here one such new test, based 
on the measurement of photon and graviton arrival times 
from the same cosmological source. As we already empha- 
sized (see ^3.1.2| l the most violent and energetic phase of the 
SMBH binary merger will be the final coalescence of the two 
BHs. Let us assume for simplicity that a luminous burst of 
light is generated at the same time as the most luminous burst 
of GWs is emitted by the coalescing binary system. If gravi- 
tons propagate at a speed lower than that of light, a delay in 
the arrival time of the burst of gravitons, relative to the burst 
of photons, may become apparent after propagation over cos- 
mological distances. 



One specific scenario in which this is exp ected is if the 
graviton possesses a non-zero rest-mass (e.g, iPauli & Fierz I 
1939).'" Assuming that Lorentz invariance is satisfied by the 
massive gravity theory under consideration, we can infer the 
bound that can be set on the graviton mass, for a given value 
of the delay between photon and graviton arrival times. Let 
us write the light travel time to a cosmological stand ard siren 
as rj (i.e. the conformal time, e.g. Eisenstein1ll997h and the 
velocity of gravitons as Cg < c, where c is the speed of light. 
By definition, we have 



c,_ Af 

C 7] 



(3) 



where Af > is the delay in graviton arrival times. The cor- 
responding Lorentz factor is 



1 



7^ 



2Af ' 



(4) 



so that from jniQC^ = hf, one could deduce the graviton rest 
mass mo from the GW frequency, /, and any measured delay, 
Af. Note that, using the conformal time ri can be expresse d 
with redshift for a given cosmology (e.g. 'Eisens teinTi 1 997h . 
and no other cosmological factors enter this calculation. 

Without any special treatment, it appears unlikely the burst 
of GWs can be isolated and thus accurately timed to better 
than about a dynamical timescale for the coalescing SMBH 
binary (i.e., one final GW cycle). This timing uncertainty will 
thus set a limit on the minimum graviton rest mass that can be 
constrained from delays in arrival times. Adopting the inverse 

Building a satisfactory Lorentz-invariant theory of massive gravity re- 
mains a challenge. Possi bly more successful Lo rentz-violating versions have 
also been proposed (e.g, IDubovskv et al.l2005l) . 
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of the frequency at ISCO as the limiting timing uncertainty, 
/isco = 2.2 X lO-'^^^'^^'Hz, for an SMBH binary with total mass 
M = 1O^''' ''M0 at z = 1, respectively, for which 77 = 1 1 .6 Gyr, 
we find that bounds on the graviton mass based on arrival time 
delays are possible as long as mo > (14,4.5, 1.4) x 10"^^ eV 
or in terms of Compton wavelength Ac < (28,89,280) pc, re- 
spectively. This bound on the graviton mass is less stringent, 
than the 1 kpc value derived from the deviations that would 
emerge in the phas ing of massive GW waveforms (iBerti et al.l 
l2005al: I Will Il2004 . In addition, any systematic delay in the 
emission of the burst of photons, which must causally follow 
the gravitational event and thus the burst of gravitons, could 
significantly reduce the quality of the constraints on mo based 
on differences in arrival times between photons and gravitons. 
In fact, in the absence of any quantitative information on such 
systematic delays at the source, one will be inevitably reduced 
to model-dependent statements, with specific emission sce- 
narios, unless observations were to reveal an unambiguous 
propagation effect (with the burst of photons significantly pre- 
ceding the burst of gravitons). 

In principle, however, a possibility exists for more accu- 
rate relative timing of the photon and graviton signals. Be- 
fore coalescence, gas present in the near environment of the 
SMBH binary would be gravitationally perturbed in such a 
way that it could radiate a variable signal with a period closely 
matching that of the leading-order quadrupolar perturbation 
induced by the coalescing binary. In this case, independently 
of the details of the electromagnetic emission process, it may 
be possible to match the periods of the electromagnetic and 
GW signals. The offset in phase between the Fourier compo- 
nents of the two signals with similar frequencies could thus 
be used to effectively calibrate the intrinsic delay in electro- 
magnetic emission at the source. Any drift in arrival-time 
with frequency between the gravitational and electromagnetic 
chirping signals, as the source spans about a decade in GW 
frequency during the last 2 weeks before merger (see Fig.|5]l, 
could then be attributed to a fundamental difference in the way 
photons and gravitons propagate on cosmological distances. 
Another immediate advantage of this frequency matching is 
that there is no need to isolate a luminous burst of GWs and 
to associate it with a corresponding burst of photons in the ob- 
served signals. Instead, the various phases of late inspiral and 
coalescence can be tracked via the GW signal, so that the rel- 
ative timing of the gravitational and electromagnetic signals 
may be known to within a fraction of the binary's dynami- 
cal time. Of course, it is difficult to assess the accuracy of this 
matching and tracking with frequency without a more detailed 
analysis. However, what is clear from our previous discussion 
is that if the delay in electromagnetic emission at the source 
can be efficiently calibrated and the absolute timing of both 
signals can be reduced to less than the binary's orbital pe- 
riod, then the direct comparison of signals could provide con- 
straints on the graviton mass, mo, of comparable or even bet- 
ter qu ahty than those obtained from the phasing of the GWs 
alone (iBerti et al.ll2Q05atlWill Il20()4 >. 

It is important to emphasize that the above limits on the 
value of the graviton mass, mo, are stron gly theory-depen dent. 
In ou r discu ssion, as in the analyses of iBerti et al.l (l2005a l) and 
IWm I (|2004|) . it was assumed that the finite value of the gravi- 
ton mass has negligible effects on the rate of decay of the bi- 
nary or other observable properties of the GW signal (such as 
polarization). Furthermore, Lorentz invariance was explicitly 
used to relate the delay in arrival times of gravitons to their 



rest mass, mo. These various assumptions can be violated in 
various theories of modified gravity, whether massive or not. 
We will mention a few possibilities of interest here in that they 
could be phenomenologically tested. 

In braneworld models with extra dimensions available to 
gravity bey ond a large, cosmological "cross-over" scale (e.g., 
iDvali et aLl l2000). from the Kaluza-Klein mechanism, mass- 
less higher-dimensional gravitons may start behaving like 4D- 
massive particles once they travel distances comparable to 
the theory cross-over scale. There would not be any single 
rest mass associated with 4D-massive gr avitons, so that the 
specific constraints discussed above or in IBerti et all (l2005ah 
would not strictly apply. That is, after cosmological propa- 
gation, the relation between arrival time and GW frequency 
may be more complex than assumed in equation (|4|i. Even 
though Lorentz invariance has been extensively tested for 
standard model fields, the possibility remains that Lorentz 
symmetry is violated in the gravity sector, especially on cos- 
mological scales. Such apparent gravity sector Lorentz viola- 
tions emerge in braneworld scenarios, if the phenomenology 
is in correctly interpr eted from a 4D-spac etime point of view 
(e.g. lChung.Kolb& Riotto 2002; Csaki et al. II2OOII) . Inde- 
pendently of higher-dimensional scenarios, specific Lorentz- 
violating theories of gravity in 4D have also been proposed 
in recent years, e.g. with a vector-tensor character (e.g. 
Bekenstein 2004; Jaco bson & M attingly 200%. 

In this context, it is thus significant that electromagnetic 
counterparts to SMBH binary mergers, together with the op- 
portunity to match and track the gravitational and electromag- 
netic signals in frequency, may offer unique tests of Lorentz 
violations in the gravity sector. It may be possible, as the 
SMBH binary decays toward final coalescence, spanning a 
range of frequencies, to measure the delays in graviton vs pho- 
ton arrival times as a function of increasing GW frequency, /. 
Therefore, the consistency of Eq. ^ with the special relativis- 
tic relation ynoc^ = h f could be tested explicitly, potentially 
revealing empirical violations of Lorentz invariance for gravi- 
tons propagated over cosmological scales. Clearly, to have 
any chance to perform such new tests of gravitational physics, 
one will need to identify the electromagnetic counterparts as 
early as possible. This may be one of the strongest motiva- 
tions behind ambitious efforts to localize these rare, transient 
events well before final coalescence. 

5. DISCUSSION AND IMPLICATIONS 

In this paper, we have considered the possibility of localiz- 
ing coalescing SMBH binaries in real time, by using the grav- 
itational waves detected by LISA to trigger searches for elec- 
tromagnetic counterparts using wide-field instruments. The 
idea presented here is that such searches may reveal a time- 
variable source - showing either quasi-periodic brightness 
variations during the inspiral stage, or a transient "burst" dur- 
ing the coalescence itself - and that this may allow identifica- 
tion of a unique counterpart. 

We first analyzed the behavior of the time-evolving local- 
ization error prior to the merger As the binary evolves in time, 
the signal-to-noise ratio, S/N, per unit time of the GW detec- 
tion increases rapidly. Therefore the measurement precision 
of physical parameters, such as the 3D position, inclination, 
and spins can improve quickly as the merger approaches. We 
have calculated the size and shape of the 3D localization ellip- 
soid, and its distribution for arbitrary binary orientations and 
sky positions, as well as how these uncertainty distributions 
evolve before the merger. We found that for typical sources at 
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z = 1-2, ~ 10 days before merger, the distance measurement 
precision reaches its Hmiting value Szwi ^ (2-10)% imposed 
by weak gravitational lensing. Redshift errors subsequently 
decouple from the sky position errors, and the 3D localiza- 
tion ellipsoid aligns with the line of sight. The sky position 
error decreases below lOdeg^, 10 days before merger, for a 
large range of masses, 10^ < M/Mq < 5 x 10^, at redshifts 
z < 1 .5. We have also found that the sky position error ellipse 
is nearly circular after the first few weeks of observation, ex- 
cept for the final ^day before the merger, especially for the 
small ellipsoids with favorable sky positions and inclinations. 
During the final day, the minor axis of the s ky ellipse shrinks 
more quickly, especially for large BH spins ( Lang & Hughes I 
l2007h . creating an eccentric sky position error ellipse. Finally, 
we found that the source position can be important: localiza- 
tion improves quickly for sources that lie in the plane defined 
by the LISA triangle, but remains much poorer in the perpen- 
dicular direction. 

We then analyzed the prospects for electromagnetic mon- 
itoring of this sky area to look for prompt counterparts and 
host galaxy candidates. It will be possible to estimate the time 
of the merger to better than an hour precision already several 
months in advance of merger. The expected number of galax- 
ies that would be present in the GW localization region is up 
to 10^ during the last weeks. However, we have shown that us- 
ing photometric redshift information and luminosity cuts con- 
sistent with the empirical black hole mass - luminosity mea- 
surements, one can decrease the number of galaxy counterpart 
candidates considerably, down to a number ^ 10^-10"^ ap- 
proximately 10 days before merger, ~ 100-1000 at 1 day, and 
^ 1-1000 at ISCO. These galaxies will be typically very faint, 
with absolute magnitudes between Mb = (-18,-16,-14) ±2 
for M = 1O'''^'^M0. Moreover, the GW localization volume 
will enclose ~ 30 quasar counterpart candidates at 10 days 
before merger, decreasing to ^ 1 quasar close to ISCO. We 
find that both the orbital inclination and the spins can be mea- 
sured with GWs to an excellent precision for each counterpart 
candidate in the field, several months before the merger. Us- 
ing this information, it may be possible to assess, in advance, 
which SMBH binaries will receive a recoil with a direction 
(and known magnitude) that is likely to plunge them into the 
surrounding gaseous disk, favoring a transient EM signal. If 
SMBH binary orbits align with the plane of the galactic disk, 
then measuring the disk orientation to ±20deg would have a 
potential for further reducing the number of counterpart can- 
didates by typically more than an order of magnitude. 

The best hope to find a merger event EM counterpart may 
be to look for variable emission produced by gas near the 
BHs. EM variability can be expected if there is ambient gas 
surrounding the coalescing binary, on the timescale of the pe- 
riodic gravitational perturbation, ~ 0.1-10 hr An electro- 
magnetic transient might take place shortly after the merger 



due to GW recoil or the sudden central rest-mass loss. The 
resulting counterpart candidates should be monitored for any 
type of EM transients or periodic variability whose period and 
time dependence is correlated with that of the known gravita- 
tional perturbation. The basic requirements of a wide field 
of view, and fast detectors, are similar to those of searches 
being planned for distant cosmological supernovae. For typ- 
ical LISA sources, a triggered EM counterpart search cam- 
paign will require monitoring a several-square degree area, 
and could aim for variability at the 24-27 mag level in optical 
bands, corresponding to 1-10% of the Eddington luminosity 
for prime LISA sources with ~ (IO^-IO^)Mq BHs at z = 1 - 2, 
on time-scales of minutes to hours. A cross-correlation of the 
phase of any variable EM signal with the quasi-periodic grav- 
itational waveform over 10-1 000 cycles may aid the detection. 
The optical emission from near the BHs may be obscured by 
large columns of gas and dust during the merger, in which 
case monitoring a similar area, to similar flux levels in X-ray, 
infrared or radio bands may be required. 

The real-time monitoring envisioned here could well make 
the difference between detecting and missing the prompt EM 
counterpart. Our results therefore also constitute a strong 
scientific case for allowing data from LISA to be optionally 
downloaded on ~ hourly timescales, to permit near-real-time 
analysis when a particularly interesting event is noticed. 

The secure identification of the EM counterpart to a GW 
event would indeed have far-reaching implications. The iden- 
tification of a counterpart would provide observational in- 
sights into poorly understood processes related to SMBH 
binary accretion and transients during a general relativistic 
merger involving the presence of gas. Another tantalizing ap- 
plication of such observations is the possibility to constrain, in 
new ways, models of gravity on cosmological scales and the 
geometry of the universe with these events. Comparing arrival 
times between the EM and GW signals may allow stringent 
bounds to be placed on the mass of the graviton and could en- 
able empirical tests of Lorentz violations in the gravity sector. 
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